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ABSTRACT 

Quiescent emission from the neutron star low-mass X-ray binary Cen X-4 is seen to 
be variable on timescales from hundreds of seconds to years, suggesting that at least 
in this object, low-level accretion is important during quiescence. Here we present 
results from recent XMM-Newton and Swift observations of Cen X-4, where the X-ray 
flux (0.5 - 10 keV) varies by a factor of 6.5 between the brightest and faintest states. 
We find a positive correlation between the X-ray flux and the simultaneous near-UV 
flux, where as there is no significant correlation between the X-ray and simultaneous 
optical (V, B) fluxes. This suggests that while the X-ray and UV emitting regions are 
somehow linked, the optical region originates elsewhere. Comparing the luminosities, 
it is plausible that the UV emission originates due to reprocessing of the X-ray flux 
by the accretion disk, with the hot inner region of the disk being a possible location 
for the UV emitting region. The X-ray/UV correlation does not favour the accretion 
stream-impact point as the source of the UV emission. 
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1 INTRODUCTION 

Neutron star low-mass X-ray binaries are often transient, 
cycling through periods of outburst and quiescence where 
the X-ray luminosity varies by over four orders of magni- 
tude between the two. In quiescence, the accretion rate is 
at less than 10 -4 of the Eddington limit. Accretion dur- 
ing quiescence should occur in a very different manner from 
the standard thin disk accretion |Shakura &: SunvaevHl973h 
during outburst. As the accretion rate drops, the inner disk 
is thought to evaporate becoming a r adiatively inefficient 
accretion flow (RIAF; see, for exam ple, iNaravan fc Yi|[l994l . 
ll995l : lBlandford fc Begelmanll 19991 ). with the accretion disk 
truncated further out. This model predicts that neutron 
stars should have a higher quiescent luminosity compared 
to black holes as the hot flow can heat the neutron star, 
where as in black holes t he energy wil l be a dvected through 
the event horizon (e.g., iGarcia et all lioOll . and references 
therein). However, the quiescent neutron star is also ex- 
pected to be hot due to heat ing of the cru st during the out- 
burst ('deep crustal heating' : I Brown et allfl99^ h and energy 
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may go into a jet rather than being advected (|Fender et al.l 
2003), so this is not a unique prediction of RIAF models. 
Thus far, there has been no strong observational evidence 
leading to a clear picture of how accretion proceeds during 
quiescence. 

One potential signature of accretion during quiescence 
is sporadic X-ray variability. Such variability in quiescent 
emission from neutron stars has been seen on all timescales 
studied - from hundreds of seconds to years t Campan a et al.l 
1991 120041 : iRutledge et all 1200 ll, [20021 : ICackett et all 12005 
20ld.[2011: Fridrikssonet all 



2011 



X-ray spectra of quies- 
cent neutron stars typically show two components - both 
a thermal, blackbody-like component and a non-thermal 
power- law component. A recent study of lon g-term quiescent 
variability in Cen X-4 jCackett et aUl201Ch confirmed that 
the thermal component, as well as the power-law component, 
is variable. Interestingly, the power-law and thermal fluxes 
appear to vary in tandem, contributing roughly the same 
fraction of the total flux at all epochs. This suggests a clear 
link between the two components during quiescence. A pos- 
sible scenario involves residual accretion onto the neutron 
star surface, where as the accretion rate onto the neutron 
star surface varies the temperature of the surface changes. 
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In order to further investigate the nature of quies- 
cent variability in neutron star low-mass X-ray binaries on 
timescales of weeks to months we observed Cen X-4 with 
XMM-Newton four times between August 2010 and January 
2011. Importantly, we also obtain optical and UV photom- 
etry with XMM-Newton to investigate the nature of broad- 
band variability during quiescence. Combining these four 
new XMM-Newton observations with one archival XMM- 
Newton observation and two Swift observations we show a 
clear correlation between X-ray and UV flux during quies- 
cence. Cen X-4 is one of very few objects where such an 
X-ray/UV study can be performed given its proximity and 
relatively low extinction. In Section [2] we describe the data 
reduction and analysis before describing our results in Sec- 
tion [3] and discussing their implications in Section [4] 



2 DATA REDUCTION AND ANALYSIS 

The four XMM-Newton observations were obtained between 
the end of August 2010 and January 2011 with each ob- 
servation lasting between 10 and 20 ksec. The ObsIDs are 
0654470201, 0654470301, 0654470401, 0654470501 and their 
respective start dates (dd/mm/yy) are 25/08/10, 04/09/10, 
24/01/11 and 31/01/11. The X-ray detectors were operated 
in full win dow mode with a m edium filter. The Optical Mon- 
itor (OM. iMason et al.|[200ll ) was operated such that data 
was collected in multiple optical/UV filters during the X-ray 
observations. In Table [1] we give details of all four observa- 
tions along with exposure times for each instrument and 
filter used. The data was reduced with the XMM-Newton 
Science Analysis Software, version 11.0.0, using the latest 
calibration files. 

In addition to the XMM-Newton data, we also analyze 
two Swift observations of Cen X-4 that also have both X- 
ray and optical/UV data. The ObsIDs 00035324001 and 
00035324002 were performed on 03/09/06 and 01/05/12. 
The X-ray detector was operated in photon counting mode. 
In the first observation, only a UVW1 exposure was taken, 
where as in the second observation exposures in V, B, U, 
UVW1, UVM2, and UVW2 were taken. Tablerjgives details 
of the Swift observations. The data was reduced using HEA- 
SOFT, version 6.11.1. In the following sections we detail the 
data reduction for the X-ray and optical/UV detectors. 

2.1 X-ray data reduction 

2.1.1 XMM-Newton 

We created calibrated event files for each observation from 
the Observation Data Files using the emproc and epproc 
commands for the MOS and PN detectors. We check for 
periods of high background by creating light curves from the 
entire detectors with 100s time binning, filtering for events 
with energies > 10 keV and PATTERN=0 for the MOS, and 
10 < E < 12 keV and PATTERN=0 for the PN. 

The background was high throughout the entire first 
observation (0654470201), with the MOS1 count rate rang- 
ing from 0.2 - 3.7 c/s, and the PN count rate from 2.6 - 
15.2 c/s. Filtering out the periods of particularly high back- 
ground would leave only a very small exposure left, thus we 
include all the available data. The other three observations 



had some short periods of background flaring at either the 
beginning or end of the observations, but the background 
was mostly low throughout. The periods of background flar- 
ing were never as bad as the first observation, and therefore 
we include all available data. Only in the first observation 
does the significant flaring there noticeably reduce the S/N 
ratio. 

For both the MOS and PN detectors, we extract the 
source spectrum from a circular region of radius 20 arcsec, 
and the background spectrum from a nearby, source-free 2 
arcmin region. The response files are generated with the ar- 
fgen and rmfgen tools, and the resulting spectrum is binned 
to a minimum of 25 counts per bin. We give net count rates 
in Table Q] 

2.1.2 Swift 

Calibrated event files were created by reprocessing the data 
using the xrtpipeline tool, and applying the standard (de- 
fault) screening criteria. We analyze data taken in photon 
counting mode. We extracted the source spectrum using xse- 
lect, and a circular extraction region with a 20 pixel radius. 
The background spectrum was extracted from an annulus 
with inner radius of 40 pixels and outer radius of 120 pixels. 
We use the xrtmkarf tool to create the ancillary response 
file, and use the appropriate response matrix based on the 
epoch of the observation and the observing mode. Given the 
shorter exposures and lower count rates compared to XMM- 
Newton we only bin spectra to 10 counts per bin. 

2.2 Optical/UV data 

2.2.1 XMM-Newton 

We observed Cen X-4 with the OM in imaging mode. How- 
ever, rather than use the standard set-up, we chose to more 
efficiently obtain exposures of Cen X-4 in 5 filters by using a 
Science User Defined mode whereby only a single window is 
used to observe just the center of the field of view (normally 
a mosaic is created from multiple windows to cover a much 
larger field of view). In this way, we were able to obtain im- 
ages in the following filters during the short X-ray exposures 
of Cen X-4 (effective wavelengths are given): V (5407A), B 
(4334A), U (3472A), UVW1 (2905A) and UVM2 (2298A). 

We reduce the OM data using the omichain tool. This 
tool processes the data with the latest calibration files, and 
then performs source detection and aperture photometry. 
The output files include both the images for each filter along 
with a combined source list containing detailed parame- 
ters for each source detected in all exposures, including the 
count rate (corrected for instrumental effects such as detec- 
tor deadtime), and AB flux densities (in erg s _1 m -2 A" 1 ). 
We visually identified Cen X-4 in all images before finding 
the corresponding source in the source list. In observation 
0654470301 Cen X-4 was not automatically detected in the 
UVM2 filter, though a faint excess in counts did seem appar- 
ent on visual inspection. We therefore ran the omdetect tool 
(part of the omichain pipeline) on its own with the signif- 
icance reduced to only requiring 1<t above the background, 
however, the source was still not detected (note the default 
is a 3<t threshold). Similarly, in observation 0654470501 the 
source was not automatically detected in UVW1 or UVM2. 
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Table 1. Details of the new XMM-Newton observations of Cen X-4, including count rates and optical/UV fluxes. AB flux densities have 



units of erg s 1 


cm 2 A. ND indicates that the 


source was not detected in 


that filter. MOS and PN net count rates are given in 


the 0.5 


- 10 keV 


energy range. The count rate is . 


given for MOS 1 only, the MOS 2 rate 


is always very 


similar. 








Detector 




0654470201 






0654470301 






0654470401 






0654470501 




or 


Exp. 


Net rate 


AB flux 


Exp. 


Net rate 


AB flux 


Exp. 


Net rate 


AB flux 


Exp. 


Net rate 


AB flux 


filter 


(ks) 


(10 _1 c/s) 


(io- 17 ) 


(ks) 


(IO" 1 c/s) 


(IO" 17 ) 


(ks) 


(10 _1 c/s) 


(io- 17 ) 


(ks) 


(10 _1 c/s) 


(IO" 17 ) 


MOS 


11.2 


0.75 ± 0.03 




19.3 


0.84 ± 0.02 




14.4 


1.20 ± 0.03 




13.4 


0.39 ± 0.02 




PN 


10.4 


3.07 ± 0.07 




16.1 


3.73 ± 0.05 




11.3 


4.21 ± 0.06 




10.5 


1.59 ± 0.04 




V 


4.0 


6.2 ± 0.6 


15.5 ± 1.4 


2.86 


9.2 ± 0.7 


23.1 ± 1.8 


1.76 


5.1 ± 0.8 


12.8 ± 2.1 


1.76 


4.8 ± 0.8 


12.0 ± 2.0 


B 


3.0 


7.0 ± 0.8 


8.7 ± 1.0 


0.0 






1.76 


8.2 ± 1.0 


10.2 ± 1.2 


1.76 


0.4 ± 1.0 


8.0 ± 1.2 


U 


1.98 


2.9 ± 0.6 


5.6 ± 1.2 


3.76 


4.6 ± 0.5 


8.8 ± 0.9 


1.76 


5.0 ± 0.6 


9.7 ± 1.2 


1.76 


2.6 ± 0.6 


5.0 ± 1.2 


UVW1 


5.0 


1.7 ± 0.2 


8.3 ± 1.1 


3.76 


1.8 ± 0.3 


8.8 ± 1.3 


1.76 


2.4 ± 0.4 


11.5 ± 2.0 


1.76 


0.42 ± 0.36 


2.0 ± 1.7 


UVM2 


5.0 


0.37 ± 0.14 


8.2 ± 3.0 


5.0 


ND 


ND 


4.2 


0.57 ± 0.15 


12.6 ± 3.3 


3.2 


ND 


ND 



Table 2. Details of the Swift observations of Cen X-4, including count rates and optical/UV fluxes. AB flux densities have units of erg 
s _1 cm~ 2 A. XRT net count rates are given in the 0.5 - 10 keV energy range. 



Detector 




00035324001 






00035324002 


or 


Exp. 


Net rate 


AB flux 


Exp. 


Net rate 


AB flux 


filter 


(ks) 


(10" 1 c/s) 


(IO" 17 ) 


(ks) 


(IO" 1 c/s) 


(10" 17 ) 


XRT 


4.5 


0.38 ±0.03 




3.8 


0.76 ± 0.05 




V 








0.30 


0.75 ± 0.07 


19.5 ±2.0 


B 








0.30 


1.11 ± 0.10 


16.5 ± 1.4 


U 








0.30 


1.35 ± 0.09 


22.5 ± 1.5 


UVW1 


6.25 


0.32 ± 0.01 


15.0 ±0.5 


0.61 


0.55 ± 0.04 


25.6 ± 1.9 


UVM2 








0.95 


0.29 ± 0.02 


24.4 ± 1.8 



Table 3. XMM-Newton obsID 0144900101, optical/UV details. 
AB flux densities have units of erg s — 1 cm -2 A. 



Filter 


Exp. 


Net rate 


AB flux 




(ks) 


(10" 1 c/s) 


(io- 17 ) 


V 


4.0 


8.3 ±0.5 


20.8 ± 1.2 


B 


4.0 


13.9 ±0.5 


17.3 ± 0.7 


U 


4.0 


6.6 ±0.3 


12.7 ±0.7 


UVW1 


12.2 


3.7 ±0.1 


18.0 ±0.5 


UVM2 


19.8 


0.81 ±0.05 


17.9 ± 1.0 


UVW2 


25.0 


0.25 ±0.07 


14.5 ± 3.8 



Once again, it did seem present (though faint) on manual in- 
spection of the image, and we therefore re-ran the omdetect 
tool with only a la threshold. This led to a detection of Cen 
X-4 in UVW1 with only a 1.9a significance, while there was 
still no positive detection in the UVM2 filter. Count rates 
and AB flux densities are given for each observation in Ta- 
ble [1] There is clearly variability seen between the different 
observations. 

We also note that only one of the two archival XMM- 
Newton observations (ObsID:0144900101) also has OM data. 
This observation took place starting on 1 March 2003, with 
exposures in V, U, B, UVW1, UVM2 and UVW2 filters. The 
data was taken in imaging mode, using the default setup 
whereby a series of 5 exposures are taken in each filter, and 
then combined together in a mosaic to fill the 17 arcmin field 
of view. During each exposure the central 2 arcmin region is 
always observed. We reduced the data in the same manner, 
using the omichain tool once again. The exposure times, 
corrected count rates and fluxes for each filter are given in 
Table H 



2.2.2 Swift 

The two Swift observati ons also include optical/UV observa- 
tions using the UVOT (|Roming et al . 2005; iBreeveld" et alj 
2010). The filters are similar, though slightly different, to 
the XMM-Newton filters. The first observation used the 
UVW1 filter exclusively, where as the second observation 
used V (5468A), B (4392A), U (3465A), UVW1 (2600A), 
UVM 2 (2246A) and U VW2(1928A) filters (wavelengths 
from iPoole et all 120081 ). Note that both the UVW1 and 
UVW2 filters are broader than the UVM2 filter, and hence 
are more sensitive. However, UVW1 and UVW2 also have 
red leaks, meaning that their sensitivity stretches to longer 
(redder) wavelengths than UVM2. We used the uvotmaghist 
tool to perform aperture photometry, using a 3 arcsec radius 
circular extraction region, and a circular, source-free nearby 
background region with a 12 arcsec radius. This tool per- 
forms photometry on all separate exposures for a given filter, 
applies the aperture correction and determines the AB flux. 
Having determined the count rates and fluxes for each indi- 
vidual exposure, we calculate the exposure- weighted average 
count rate and flux, as given in Table [2] 



2.3 X-ray Spectral analysis 

We fit the X -ray spectra using XSPEC version 12.7.0 
i|Arnaudlll996l ). Each observation is fit separately, but when 
fitting the individual XMM-Newton observations we jointly 
fit the MOS 1, MOS 2, and PN spectra, with the param- 
eters the same for all detectors. We fit the spectra using 
a neutron star atmosphere plus a power-law, all modified 
by Galactic photoelectric absorption. We use the nsatmo s 
model for the neutron star atmosphere l|Heinke et al.ll2006h . 
This is the same model we used to fit six previous obser- 
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vations of Cen X-4 (jCackett et al.ll2010h . thus allowing for 
a direct comparison with these previous result^]. For the 
same reason, we choose to fix Nh — 4.9 x 10 20 cm -2 , the 
va lue we found fr om jointly fitting six previous observations 
in lCackett et al.l (|2010h . Also, to allow a direct comparison 
with our previous re sults, we fix the distan ce to 1.2 kpc 
Chevalier et al.lll989T ) in the spectral fits. In lCackett et al.l 
20101 ) we investigated which parameters were variable, find- 
ing that the thermal component and the power-law compo- 
nent must both vary. For the thermal component, either the 
effective temperature, emitting radius or both can vary, and 
we found that equally good fits are achieved regardless of 
which parameter is left variable. Here, for simplicity, we fix 
the neutron star radius and allow the temperature to vary 
between epochs. This choice does not affect the results since 
the thermal component can be fit equally well with either a 
variable temperature or a variable radius, and it is only the 
flux that we are mostly concerned with here. We assume a 
neutron star with R = 10 km and M — 1.4 Mq, and that the 
entire surface is emitting. We fit the spectra in the 0.5 - 10 
keV range and best fitting parameters are given in Table [4] 
Uncertainties are quoted at the la level throughout. 




Years since 1 980 

Figure 1. Long-term X-ray quiescent light curve of Cen X-4, 
showing variability on all timcscales. (a) Unabsorbed 0.5 - 10 
keV flux, (b) Effective temperature for an observer at infinity. 



3 RESULTS 

The long-term X-ray quiescent light curve of Cen X-4 is 
shown in Figu re [Q where we include the observations from 
ICackett et all HqIO) as well as the four new XMM-Newton 
and two Swift observations presented here. Notably, signif- 
icant variability is seen on the weeks to months timescales 
probed by these new observations. Most significant is the 
factor of 3 drop in flux between the third and forth obser- 
vations taken at the end of January 2011 and separated by 
only 7 days. 

It is clear from looking at the OM count rates and fluxes 
(Table [T]) that Cen X-4 is also highly variable in the op- 
tical/UV part of the spectrum during quiescence. In Fig- 
ure [5] we plot the optical/UV fluxes against the 0.5 - 10 
keV unabsorbed X-ray flux. In order to compare the ener- 
getics, we convert the flux densities (erg s _1 m -2 A -1 ) that 
are the output from the omichain to fluxes by multiplying 
by the FWHM of the filter bandpass. We use: 684 A (V), 
976A (B), 810A (U), 620A (UVW1), 439A (UVM2) for the 
XMM-Newton filters, and 769 A (V), 975 A (B), 785 A (U), 
693A (UVW1), 498A (UVM2) for the Swift filters. We also 
de-redd en the UV/optical fluxes. W e use the gas-to-dust ra- 
tio from iPredehl fc Schmittl (|l995l ) in order to convert from 
equivalent hydrogen column to E(B — V). We then use the 
interstellar extinction curve of lCardelli etHl (|l989l ) to es- 
timate the extinction correction at the wavelengths of each 
filter. 

We test for a correlation between the optical/UV 
and X-ray fluxes using a simple Pearson's linear corre- 
lation test. The correlation coefficients we find are r = 
0.553,0.778,0.982,0.974,0.977 for the V, B, U, UVW1 and 
UVM2 filters respectively. Given the number of data points, 
and a two-tail test (no a priori knowledge of positive or neg- 
ative correlation), this corresponds to a positive correlation 

This previous analysis looked at data from a range of different 
missions (1 ASCA, 2 Chandra, 2 XMM and 1 Suzaku) 



at the 0.745, 0.879, 0.99952, 0.99979, 0.977 confidence levels. 
Thus, the U and UVW1 fluxes are correlated with the X-ray 
flux at greater than 3<r, while the UVM2 flux is correlated 
at greater than 2o. There is no significant correlation be- 
tween the V or B fluxes and the X-ray flux, indicating that 
this is a UV only correlation. We also look to see if there 
is a correlation between UVW1 and V band fluxes, finding 
r = 0.63, which corresponds to a positive correlation at the 
0.816 confidence level. Thus, there is no significant correla- 
tion between the V band and UVW1 fluxes. 

We now concentrate on the UVW1 versus X-ray cor- 
relation further, given that it is the most significant cor- 
relation. Obviously, given the high linear correlation coef- 
ficient, a simple straight-line fits the UVW1 - X-ray cor- 
relation well. We get the following best-fitting parameters: 
fuv = (0.16 ± 0.01)/x - (0.57 ± 0.20), giving a reduced- 
X 2 = 1.29. We also test fitting a power-law of the form 
fuv = a(fx — fxo) b , where fxo is a constant to allow a non- 
zero X-ray flux when fuv ~ 0.0, or vice versa. The best-fit 
is very close to a linear relationship, with b — 1.02+Q14, and 
giving a reduced-x 2 = 1.61. With only 7 data points we can 
clearly not make any strong conclusions about the form of 
the correlation, especially as there is only one observation 
at a high luminosity, where a deviation from a linear rela- 
tionship would become apparent. Further data is needed to 
strengthen and define the correlation. 

We also look for whether there is an X-ray/UV correla- 
tion on shorter timescales by looking at the UV and X-ray 
light curve during the two Swift observations. The is possi- 
ble with Swift because the observations consist of multiple 
~ 1 ksec exposures, each giving separate UV images. With 
XMM-Newton however, there is only one exposures in each 
filter during the entire observation and so we cannot look 
for X-ray-UV correlations during the XMM-Newton obser- 
vations. The first Swift observation consists of 6 Swift or- 
bits of data, and the UVW1 filter was used for all of each 
orbit. We can therefore compare the X-ray and UV light 
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Table 4. Spectral fit parameters. The flux is given in units of erg s — 1 cm -2 . The column density was fixed at TVh = 4.9 X 10 20 cm -2 
in all spectral fits. We also assumed a neutron star radius of iO km, and mass of 1.4 Mq, and a distance to Cen X-4 of 1.2 kpc. The 
power-law normalization is defined as photons keV - 1 cm~ 2 s _1 at 1 keV. We define the thermal fraction is the ratio of the unabsorbed 
0.5 - 10 keV thermal flux to the total unabsorbed 0.5 - 10 keV flux. 



XMM-Newton Swift 



Parameter 


0654470201 


0654470301 


0654470401 


0654470501 


00035324001 


00035324002 


fcT<g (eV) 


55.9 ±0.7 


56.5 ±0.3 


60.8 ±0.3 


49.0 ±0.6 


62.4 ±2.3 


65.2 ±5.7 


Power-law index, T 


1.77 ±0.21 


1.62 ±0.10 


1.38 ± 0.10 


1.94 ±0.19 


1.51 ±0.86 


1.79 ±0.53 


Power- law norm (10 — 5 ) 


6.6 ±1.4 


6.3 ±0.7 


5.5 ± 0.7 


3.6 ±0.7 


5-8_ 4 ; 8 




Unabs. 0.5 - 10 keV flux (10~ 12 ) 


0.92 ±0.08 


1.00 ±0.02 


1.31 ± 0.03 


0.44 ±0.01 


1.38 ±0.38 


2.85 ±0.30 


Thermal fraction (%) 


58 


56 


62 


59 


68 


41 


xl (dof) 


0.91 (106) 


0.93 (168) 


0.99 (129) 


1.27 (71) 


1.31 (13) 


0.67 (24) 



-I 1 1 1 1 1 1 1 1 1 1 1 1 1 1— 



-I 1 1 1 1 1 1 1 1 1 1 1 1— 




10 20 30 

0.5 - 10 keV unabsorbed flux (10~ 13 erg s" 1 cm" 2 ) 

Figure 2. Optical and UV fluxes versus 0.5 - 10 keV unabsorbed 
flux from XMM-Newton (black ) and Swift (red). A clear corre- 
lation between the X-ray and the U, UVW1 and UVM2 fluxes is 
apparent. The solid line shows the best-fitting straight line to the 
UVW1 /X-ray points. 



£ 0.1 



0.3 g 



1.5x10* 2x10* 
Time (s) 



2.5x10 4 3x10* 



ft 



Time (s) 

Figure 3. Swift X-ray (black) and UVW1 (red, triangles) light 
curves during the first (top) and second (bottom) Swift obser- 
vations (00035324001 and 00035324002). The XRT data has 200 
second time binning. 



curves on timescales of a few thousand seconds. The sec- 
ond Swift observation comprised three orbits, with a UVW1 
exposure during each orbit. For both observations, we cre- 
ate a background-subtracted Swift/XKT light curves, using 
the xrtlccorr tool to perform count rates corrections for vi- 
gnetting, PSF etc. The UVW1 count rates are determined 
using the uvotmaghist tool (as described above). We show 
the light curves in Figure [3] The XRT and UVW1 count 
rates are generally well correlated during both observations. 

The XMM-Newton observations only perform one ex- 
posure in each filter during a given observation. Therefore, 
while the X-ray flux we obtain is from the entire observation, 
the optical/UV fluxes are only obtained during some fraction 
of the observation. The X-ray exposures range from approx- 
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imately 10 - 20 ksec, where as the optical/UV exposures 
are between 1.76 and 5 ksec. Variability during the observa- 
tion could therefore add scatter to the UV and X-ray corre- 
lations. Examining the X-ray lightcurves from the XMM- 
Newton observations there is variability present (as has 
also been seen during previous XMM-Newton and Chandra 
obser vations, see, e.g., ICampana et al ] |2004l : ICackett et al.l 
2010). However, any subset of 2 ksec or longer will not de- 
viate too far from the mean X-ray flux, and thus not add a 
large amount of scatter. 



4 DISCUSSION AND CONCLUSIONS 

Our observations of Cen X-4 during quiescence have shown 
that the X-ray and near-UV emission are correlated, where 
as there is no significant correlation between the X-ray and 
optical emission. The X-ray/UV correlation is seen both on 
timescales of weeks - years, as well as during observations 
spanning less than one day. In order to investigate the phys- 
ical origin of this correlation, it is important to consider the 
energetics - if the correlation is due to reprocessing of the X- 
ray emission, then the near-UV flux must be only a fraction 
of the X-ray flux. 

We find that the UVW1 flux is approximately 10% of 
the 0.5 - 10 keV flux. However, much of the X-ray flux from 
the hot neutron star surface is emitted below 0.5 keV, and 
thus it is important to take this flux into consideration. 
Furthermore, the power-law component will likely extend 
beyond 10 keV, though its shape there is uncertain. To pro- 
vide an estimate for the total X-ray flux, we extrapolate the 
best-fitting neutron star atmosphere plus power-law models 
to cover the range 0.01 - 100 keV. For the faintest obser- 
vation (XMM-Newton: 0654470501) we find a factor of 3 
increase in the flux, where as for the brightest observation 
(Swift: 00035324002) we find a factor of 2.6 increase. The 
higher temperature during the brighter observation means 
more of the neutron star atmosphere flux is in the 0.5 - 10 
keV band, hence the slightly lower increase in flux when ex- 
tending the energy range. Of course, the near-UV emission is 
in more than just the UVW1 filter. Combining the U, UVW1 
and UVM2 fluxes, we get an increase by a factor of 4.0 and 
2.7 for the faintest and brightest observations respectively, 
over using just the UVW1 flux. Thus, when accounting for 
emission over a wider wavelength range, the near-UV flux is 
still approximately 10% of the X-ray flux. Such a fraction of 
the X-ray flux is reasonable for r eproce ssed emission. 

van Paradiis fc McClintockl (|l994f ) consider X-ray re- 
processing in the accretion disk and show that the repro- 

1/2 

cessed emission should be proportional to L-^ R, where R 
is the outer radius of the accretion disk. We, however, find 
that the UV flux increases approximately linearly with the 
X-ray flux, though, with the limited number of observations 
and range in flux the exact relationship is uncertain. Con- 
sidering the relevant timescales for reprocessing, the light- 
travel time from the neutron star to the outer disk is less 
than one second, and so we would not expect to see any time 
delay between the X-ray and UV lightcurves. This is what 
we observe (Fig. [3j where the X-ray and UV vary in tan- 
dem during Swift observations. The viscous timescale from 
the outer disk in these objects is significantly longer (of the 



order of weeks), and therefore cannot be used to explain the 
simultaneous variability in UV and X-rays. 

Another source where a similar correlation has been 
found is the quiescent black hole V404 Cyg, where a corre- 
lation bet ween X-ray and H a emission has been previously 
observed ()Hvnes et al.ll2004h . with a weaker correlation be- 
tween the X-ray flux and the optical continuum. Irradia- 
tion/reprocessing can readily explain that correlation as the 
entire Ha line is seen to respond quickly to X-ray variabil- 
ity. Here too, we see that the X-ray and UV flux seem to 
respond simultaneously during the Swift observations (see 
Fig. [3]), though we do not have the data to be able to es- 
tablish any lag. This prompt response also supports irradi- 
ation/reprocessing as the origin of the correlation observed 
here. However, we do not observe any clear correlation be- 
tween the optical continuum and the X-ray flux (unlike in 
V404 Cyg), which would suggest that the UV and optical 
emitting regions are well-separated. 

The optical emission should be dominated by the donor 
star, as well as some contribution from the cooler outer re- 
gions of the accretion disk. The UV emission is too hot 
to arise from the donor star, and the X-ray/UV correla- 
tion suggests it must arise from a region that can see the 
central X-ray source, such as the inner region of the ac- 
cretion disk. The correlation between X-ray and UV flux 
would seem to rule-out the accretion stream-impact point 
as the location for t he UV emission, which had be en pre- 
viously suggested bv lMcClintock fc Remillardl (2000) based 
on the Mg II A2800 emission line in an HST/STIS spectrum 
of Cen X-4 not showing a double-peak as would be e xpected 
for emission from the disk. iHvnes fc Robinson! (|2012l ) closely 
examine the UV SEDs of 3 quiescent black hole and 1 quies- 
cent neutron star and also rule out the stream-impact point 
as the origin of the UV emission. They find that the mass ac- 
cretion rate required to give the observed UV emission at the 
stream-impact point in the black hole GU Mus and neutron 
star Aql X-l is 10 times higher than realistic average rates 
from these source based on their accretion histories. They 
therefore favour an origin for the UV emission as being lo- 
cat ed in the hot inner regio n of the disk as also suggeste d 
bv lCampana fc Stellal (|2000h and iMcClintock et alj (|2003h . 
Our results also support such a picture. 

It is important to consider whether there is a viable 
alternative to X-ray reprocessing as the origin of the X- 
ray/UV correlation. For instance, could the emission in 
Cen X-4 be jet-dominated and hence the UV flux and X-ray 
power-law both originate from the jet? While radio emission 
has never been detected from a quiescent neutron star, it has 
been detected in quiescence in black hole X-ray binaries. For 
instance, a ra dio detection of th e black hole hole A0620— 00 
in quiescence (jGallo et al . 2006) implies that there is a radio- 
emitting outflow during quiescence. Broadband modeling of 
the quiescent SED in A0620— 00 shows that the emission can 
be fitted with a maximally jet-dominated model with the 
jet e mission dominatin g from radio through to the soft X- 
ravs dGallo et al]|2007h and implies there are strong outflows 
IjFroning et al 1l201ll ). Similarly, V404 Cyg has been detected 
at radio wavelen g ths in quiescence dHiellming et al.l |2000| ; 
iGallo et al.1 120051 ; IHvnes et af]|2009h . with the radio spec- 
tru m consistent with synchrotron emission from a steady 
jet (|Gallo et al.ll2005h . In fact, it has been shown that quies- 
cent black hole X-ray binaries could all be in jet-dominated 
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states ijFender et al.l l2003h . Neutron stars X-ray binaries, 
on the other hand, are known to be a factor of 30 or so 
fainter radio sources than black hole binaries durin g out- 
burst (|Fender fe Kuulkersl l200ll : iMigliari et all l2003h . and 
thus, if they follow the same dependence between the X-ray 
luminosity, Lx, and radio luminosity, La, then neutron stars 
could als o be jet-dominated although at lower Eddington 
fractions (|Fender et al.l [20031 ). However, neutron stars have 
been shown to have a different dependence between Lx and 
Lr, implying that they neve r reach a jet-dominated state 
and remain X-ray dominated (|Migliari fc Fenderl200fJ V This 
therefore suggests that jet emission is unlikely to account for 
the X-ray/UV correlation observed here in Cen X-4. 

In conclusion, we have observed a significant X-ray/UV 
correlation in Cen X-4, whose most likely explanation is due 
to irradiation of the hot inner accretion disk by a central 
X-ray source leading to reprocessed UV emission. Further 
X-ray/UV monitoring will help more firmly establish the 
link between X-ray and UV emission in Cen X-4. 
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